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Abstract. The direct linearization structure is presented of a "mild" but significant general- 
ization of the lattice BSQ system. Some of the equations in this system were recently discovered 
in [J. Hietarinta, J. Phys A: Math. Theor. 44 (2011) 165204] through a search of a class of 
three-component systems obeying the property of multidimensional consistency. We show that 
all the novel equations arising in this class follow from one and the same underlying structure. 
Lax pairs for these systems are derived and explicit expressions for the A'^-soliton solutions are 
obtained from the given structure. 



1. Introduction 
The lattice Boussinesq (BSQ) equation 



p' - q-' 



P — q + Un+l,m+l — 'Un+2,m P — Q + Un,m+2 — "Un+l.m+l 

= {p-q + Un+l,m+2 - Un+2,m+l){'^P + 9 + Un,m+1 " Un+2,m+2) 

-{p-q + ){2p + q + ) , (1-1) 

appeared in |24] as a first higher-rank case in what was coined the "lattice Gel'fand-Dikii (GD)" 
hierarchy of equations. These are equations on the space-time lattice, labelled by an integer N 
associated with the order (or equivalently the number of effective components) of the system. 
Whilst for N = 2 the equation in this class is the lattice potential Korteweg-de Vries (KdV) 
equation, which is an equation associated with a quadrilateral stencil of vertices, the lattice BSQ 
appears in this class for = 3 and is defined on the following 9-point stencil: 



u 



u 
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The notation employed in (II. ip is illustrated in this Figure: u = Un,m denotes the dependent 
variable of the lattice points labelled by {n,m) £ 1?; p, q are continuous lattice parameters 
associated with the grid size in the directions of the lattice given by the independent variables 
n and m respectively, and for the sake of clarity we prefer to use a notation with elementary 
lattice shifts denoted by 

Thus, as a consequence of this notation, we have also 

U = Un+l,m+l ■ U = Un+2,m+l , U = Un+l,m+2 , U = Un+2,m+2 ■ 

Together with the lattice BSQ, also lattice versions of the modified BSQ (MBSQ) and Schwarzian 
ESQ equations have been found, cf. |24lll81[T9t l30| as well as their semi-continuum limits (leading 
to differential-difference analogues of these equations). The lattice BSQ equation has (re)gained 
considerable attention in recent years, w.r.t. similarity reductions \30\ I28j . soliton solutions 
[T2I [T7] and its Lagrangian structure [16] . 

It was shown in [30], cf. also [28], that eq. (jl.ip can also be written as the following coupled 
system of 3-field equations (for fields u, no,i and ui^q) around an elementary quadrilateral: 

uo,i + uifl = puQ - puQ + uqUq , (1.2a) 
uo,i + uifi = quo - QUO + uqUo , (l-2b) 

— — =p'^ +pq + q^ - {p + q){uo - uo) - uouo + Si,o + "Uo,! • (l-2c) 

p-q + UQ-uo 

In the recent paper [IT] a one-parameter generalization was presented of this system (after a 
point transformation on the dependent variables), as well as of two other 3-component BSQ 
systems, comprising the lattice modified BSQ (MBSQ) and Schwarzian BSQ (SBSQ) systems, 
through systematic search based on the multidimensional consistency property. Subsequently, 
in |13j soliton solutions for this extended BSQ system were derived. The main structural feature 
that emerged from the latter study was that the dispersion relation defining the relevant discrete 
exponential functions involve a cubic of the form 

uj^ -k^ - P{uj^ - A;2) = , 

where k is like the wave number of the solution, whilst uj the corresponding angular frequency. 
Noting that the parameter /3 corresponds to the deformation of the 3-component system, the 
limit /3 — 7> would lead to the original lattice BSQ system, in which case the dispersion would 
involve a cube root of unity for the parameter co (i.e., uj exp(27ri/3)A; in this limit). 

In the present study, we generalize the latter relation to an arbitrary cubic, i.e. we consider 
a dispersion of the form 

G{uj, k) := uT' + a2U? + aiw - {k^ + a2k'^ + aik) = , (1.3) 

and show that not only the extended lattice BSQ of [11], but also the generalizations of the 3- 
component systems of lattice MBSQ and SBSQ, which remained as yet unidentified, arise from 
one and the same structure. We will derive this structure, in the case of the extended dispersion 
relation (II. 3p . following the direct linearization scheme developed for the case that ai = 02 = 
in |24] . and subsequently derive from this the main equations as well as the relations between 
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the various fields arising in the extended lattice BSQ systems of [H]. In particular we derive 
the following coupled system 

VaWb = ^ p~ p ^ ^ , (l-4a) 

VaWb = PaU - PbU , VaWb = QaU - QbU , (l-4b) 

in which the quadrilateral multilinear function Q is given by 

Qp^q{u,U,U,u) \= PaPb{uU + Uu) - QaQb{uU + U^^ , (l-4c) 

and in which 



Pa = V-Gi-P, -a) , Qa = ^-G{-q,-a) , G{p,q) = - + a2{p^ - q^) + ai{p - q) . 
The system, by elimination of Va and Wb leads to the following 9-point equation for u 

%,q{u, n, U, u) {QaU - QbU) (PaU - PbU) {QaPbU - PaQbU) 

— " z — s z ■ ^ 

Qp,giu, U, U, U) (p^^ _ {QaU - QbU) {QaPbU - PaQbU) 

Eq. ()1.5p is the extended form of the BSQ analogue of the (Q3)o equation, and which reduces 
to the analogous equation given in [21] by setting ai = Q2 = 0. It is related by a point 
transformation to an equation given already in [191 [30] , and which constitutes the rank 3 analogue 
of the NQC equation of [25] (associated with the lattice KdV equations). In [21] the derivation 
of a BSQ analogue of the full Q3 equation in the ABS list was presented, but it gives rise to a 
rather complicated coupled system. Thus, even though those results can in principle be readily 
extended to the more general dispersion relation (|1.3p as well, we will refrain from doing so in 
the present paper. As a concrete spin-off of the derivations presented here, we derive also the 
Lax pairs for the extended BSQ system, whilst A^-soliton solutions in the Cauchy matrix form 
can be readily inferred from the direct linearization structure. 

2. Direct linearization and constitutive system 

In this section we present the constitutive relations of the direct linearization scheme for a 
class of lattice systems generalizing the lattice Gel'fand-Dikii (GD) hierarchy. The members of 
this hierarchy are labelled by an integer iV, which in the case of [24] were associated with the 
corresponding roots of unity uj = exp{27ri/N), the case of iV = 2 corresponding to the lattice 
KdV systems, whilst the case = 3 yielding the lattice BSQ systems. In the present paper the 
root of unity is extended to a dispersion function ujj{k) associated with the roots of arbitrary 
fixed A^**^ order polynomials. 

2.1. Infinite matrix scheme. The starting point of the construction is the following integral 

N 

C = X]/ , (2-1) 

.7=1 -^^J 
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over a yet unspecified set of contours or arcs Tj in the complex plane of a variable k, and in 
which for > 2 a positive integer, ujj{k){j = 1,2, •• • , — 1) and ujNik) = k are the distinct 
roots of the following algebraic relation 

N 

G{u}, k) := g{uj) - g{k) = , where g{k) = ^OjA;^' (2.2) 

i=i 

is a monic (i.e. oat = 1) polynomial with coefficients aj. The factors p^, a^' are discrete 
exponential functions of k, k' respectively, given by 

Pfc(n,m) = (p+fe)"(g + A;)™pfc(0,0) , au'{n,m) = {p - k')-^{q - k')-'^ak,{<d,Q) , (2.3) 

the infinite-component vector = {k^)j^z denotes the column vector of a basis of monomials 
in the variable k, and c^, denotes the transpose of c^/. The integration measures dXj{k) remain 
unspecified, but we assume that basic operations (such as differentiations w.r.t. parameters or 
applying shifts in the variables n and m) commute with the integrations. 

In accordance with the notation introduced in section 1, we denote shifts over one unit in the 
variables n and m respectively by the operations ~ and ^, which implies for the factors pk, (Tj.' 

Pk = Pk{n+l,m) = {p + k)pk , pk = Pkin,m + l) = {q + k)pk , 

and 

ak' = crk'{n + l,m) = {p- k'y^ak' , CTfc' = o-fe'(n, m + 1) = (g - A;')~^o-fc' • 
These relations imply the following linear relations for the matrix C 

C{p-'K) = {p + K)C , C{q-'K) = {q + K)C , (2.4) 

where we have introduced the matrices A and *A which are defined by their actions on the 
vector Cfc, and on its transposed vector, as follows: 



Besides, from the definition of C together with (12. 5p . one has 

g{A)C = Cg{-'A) , 
and as a consequence, using the fact that g{uj) — g{k) = Y\^^i{oj — ojj{k)) , also 



N 



J[{u:,{-p) - A) 

N-l 

nK(-p)-A) 



N 



(2.5) 
(2.6a) 

(2.6b) 



C = C\[{u:^{-p)+'K) , 

leading to 

N-l 

C = C \[{u:,{-p)+'K) , (2.6c) 

(and a similar relation to (]2.6cP with p replaced by q and C by C). The following ingredients 
determine the structure: 

i) A. Cauchy kernel Q, defined by the relation fiA + *Ar2 = O , where O is a rank 1 projection 
matrix, obeying = O. 

a) A. matrix U obeying the relation U = C — U flC. 
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in) A vector Uk defined by Uk = Pk{ck — UflCk)- 

In terms of these objects the fohowing sets of relations involving the shift ~ can be derived: 



U 



U{p- *A) 

'n-1 



{p + A)U -UOU , 

N-l 



(2.7a) 
(2.7b) 



N-2 



j=0 



n(u;K-p)+*A) 



.1=1 



O 



N-l 



n (^/(-p)-A) 



l=j+2 



u 



N 



N 



N 



i-1 



oijA^ U -U^aj Yjy- *A)'OA-''-i-' U 



(2.7c) 



The derivation of ()2.7p follows the same procedure as in [21], and they reduce to the relations 
obtained there in the special case that ujj{k) = exp{27ri j /N)k, i.e. the case that all coefficients 
aj = 0, j = 1, . . . , N — 1 in ()2.2p . By virtue of the covariance of the dynamics in terms of the 
variables n and m, similar relations to (j2.7ap and ()2.7bp hold with the shift ~ replaced by ^ 
while replacing p hy q. The associated linear problems (Lax pairs) are derived in terms of the 
object Uk, for which we have the following set of constitutive relation, [24j . 



N 



Ylico^i-p) - k) 



Uk 



Uk 



{p + A)uk -UOuk , 



N-l 



n(^,(-p)-A) 



Uk 



(2.8a) 
(2.8b) 



N-2 



.1=1 



N 

y^^ajk^Uk 



i=o 

N N i-1 



O 



N-l 



n (^/(-p)-A) 



l=j+2 



Uk , 



ajA^ Uk-UY «i *A)'OA^-i"' Uk 

j=l j=l 1=0 



(2.8c) 



and similar relations to (j2.8ap and ()2.8bp with the shift replaced by ^ while replacing p by q. 
The relation ()2.7bp can be easily derived from the formal relation 



U 



g{- p)-9{-'A) 
—p + *A 



5-(A) - g{-p) 



p + A 



U 



+u 



^ g{-p) - g(A) _ g(-p) -g(-*A) ^ 



P + A 



p 



whilst (|2.7cp follows from 

Ugi- *A) = giA)U - Ui^giA) - g{- 'A)n)U . 
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Similarly, the linear relation (I2.8b[) for follows from 

g{A) -g{-p). 



G{k, -p)uk 



p + A 



-Uk 



U 



ff(-p) -g(A) _ g(-p) -g(-*A) ^ 



p + A 



P 



whereas (|2.8cp follows from 

gik)uk = g{A)uk - U{ng{A) - g{-*A)n)uk . 

For any fixed value of A^, these abstract equations form an infinite set of recurrence relations 
defining the dynamics in terms of the independent variables n, m on the objects U and Uk 
taking values in an abstract vector space V, and an adjoint vector *e in its dual V*. To give a 
concrete realization, we can choose a fixed vector e and use the matrices A and *A to define a 
natural gradation in this vector space. In terms these, we realize quantities such as C and U as 
infinite-dimensional matrices. In fact, setting 



*eA*C/*A-'e 



(2.9) 



we obtain infinite, i.e., 'L'x'L matrices, on which A and *A act as index-raising operators acting 
from the left and from the right. In this realization the matrix O = e *e is the projector defined 
by these vectors, and through the relations 



e Cfc = 1 



1 



(2.10) 



this defines the "central" entry in the space of Z x Z matrices. We will define concrete objects 
in terms of the matrix U , where we allude to this realization. 



3. Lattice BSQ case 

In this section we will now apply the structure exhibited in section 2 to the case = 3, which 
in |24j is the one leading to the BSQ class of systems. We recall that the case N = 2 leads to 
the KdV class of lattice systems, but in that case it can be shown that the extended dispersion 
relation does not lead to new results. However, in the BSQ case {N = 3) it does. We restrict 
ourselves in this paper primarily to this case of = 3, but it is clear that the general system 
derived in section 2 can be studied in fairly similar ways for any (integer) N > 2. Most of the 
relations given in this section reduce to the ones given in [30] setting ai = a2 = 0. 

3.1. Basic objects and their relations: BSQ case {N = 3). Let us define 

G{lo, k) := Lo'^ -k^ + a2(u? - k^) + ai(w - k). (3.1) 
Obviously, for = 3 the algebraic relation ()2.2p reduces to 

G(a;,A;) = 0, (3.2) 
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and its roots are denoted by uji{k), uj2{k) and u}3{k) = k. The set of relations ()2.7p take the 
form: 



t/(p-*A) = {p + A)U-UOU, 
U (uii-p) + 'A)iu2i-p) + *A) = (uii-p) - A){u;2i-P) - A)U 

+U [{p - a2)0 - (O A - 'AO)] U , 

3 3 3 j-1 

U aj{- 'Ay = U-UY,(^j Z](" *A)'OA^-i-' U ,(3.3c) 



(3.3a) 
(3.3b) 



and similar relations to ()3.3ap and ()3.3bp with the shift replaced by ^ while replacing p by q, 
as well as 



Yiiu^{-p)-k) 



Uk = {p + A)uk -UOuk , (3.4a) 
Uk = (ujii-p) - A){u)2{-p) - A)uk 

+U[{p-a2)0-{0A-'A0)]uk, (3.4b) 

3 3 3 j-l 

Y ajk^Uk = Y Uk-UY «i ^(- *A)'OA^-i-' Uk , (3.4c) 

j=l j=l j=l 1=0 

(and similar relations to (j3.4ap and (j3.4bp with the shift ~ replaced by ^ while replacing p by 
q). For the sake of obtaining from the constitutive relations (j3.3p closed-form equations, we 
introduce the following objects: 



Va 
Sa 



l-'e{a + A)"^Ue , Wb := 1 + 'eU {b - 'A)'^ e , 

a- 'e{a + A}-^U 'Ae , t,, := -6 + *e A J7 (6 - *A)~^ e , 



,2 i. 



'e (a + A)-^ U'A^e , := + 'e A' U {b - 'A)'' e 



2 , t„ A 2 



(3.5a) 
(3.5b) 
(3.5c) 



but in particular the following object: 



Sa,b ■■= *e (a + A)-i U {-b + 'A)-^e. 



(3.6) 
(3.7a) 



Then for these objects and Uij defined in (j2.9p we have the following relations 

{p^ - a2P + ai)uij + {p- a2)uij+i + Uij+2 

= {P^ - OL2P + ai)MiJ -{p- a2)Mi+l,j + Ui+2j + {P- OL2)UifiU0j - UifiUl^j + Ui^lUQj , 

(3.7b) 

as well as 

l + {p- a)Sa,b - {P- b)Sa,b = VaWb , (3.8a) 
{p + a + b-a2)+ PbSa,b - PaSa,b = SaWb - Vah + (p - a2)VaWb , (3.8b) 
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and 



and 



and 



Sa = (p + Uo)Va - (p - a)va , (3.9a) 

^6 = {p- b)wb - ip- uo)wb , (3.9b) 

ra = PSa - (P- a)Sa + VaUQ,! , (3.10a) 

Zb = {p - b)tb - ptb + uifiWt , (3.10b) 



ra = PaVa - (P - UO - a2)Sa- Hp - 0!2)ip - Uo) +Uifi + ai)Va , (3.11a) 

Zb = PbWb + {p + uq - - {{p - 02) {p + uq) + uo^i + ai)wb , (3.11b) 

where we have used a shnphfied notation uq = no,o and Pa, Pb are defined as 

G{-p,-a) . 2 2n / N G{-p,-b) , 2 , ,2n / ,\ 

Pa = = (P +ap+a )-a2ip+a)+ai , pb = — ; = ip +bp+b )-a2{p+b)+ai . 

a — p b — p 

(3.12) 

All relations (j3.8p - (j3.11|) . which are in a slightly different form have already appeared in |30)n . 
also hold for their hat-g counterparts obtained by replacing the shift ~ by the shift ^ whilst 
replacing the parameter phy q. 

From equations (j3.7ap . (j3.9ap . (|3.11ap and their hat-g counterparts one can get a further set 
of relations 

{p + q-Uo + — - a2){p - q + Uf) - Uq) = Pa— - qa— , (3.13a) 

Va Va Va 

p - q + uo -uo = {p - a)^ - {q - a)^ , (3.13b) 

Va Va 

and similarly, from (j3.7ap . ()3.9bp and (|3.11bp . 

[p + q + UQ- — a2){p - q + Uq - Uq) = pb-^ qb— , (3.14a) 

Wb Wb Wb 

p - q + uo-UQ = {p -b)— - {q-b)— . (3.14b) 

Wb Wb 

3.2. Closed-form lattice equations. One can immediately get two closed-form lattice equa- 
tions. One is composed of (|3.9ap . its hat-g counterpart and (|3.13ap . i.e., 

= (p + uofva - (p- a)va , Sa = (q + Uo)Va - (q - a)va , (3.15a) 
{p + q -uo + — - a2){p - q + uq - Uq) = Pa— - qa— ■ (3.15b) 



^We note in passing that the yet unidentified case (A) of 3-component BSQ type systems which was presented 
in the recent paper [11] can be identified with eqs. (5.3.7a) together with (5.3.14) of [30], where, up to a point 
transformation, x can be identified with Va, z with u and y with Sa- By duality this case of [TT] can also be 
identified with eqs. (5.3.7b) and (5.3.16) of [30] identifying x with wp, z with u and y with tp. 
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Another is composed of (I3.9bp . its hat-q counterpart and (]3.14b[) . i.e., 

tb = (p- b)wb - (p- uo)wb , tb = {q- b)wb - {q - uo)wb , (3.16a) 

{p + q + uo - ^ - a2){p - q + uo -uo) = pb^ - qb^ ■ (3.16b) 

Wb Wb Wb 

To get further closed form lattice systems we take i = j = in (j3.7p and its hat-g counterpart 
and get 

puQ - Mo,i = PUQ + uifl - uqUq , quo - uq,! = quQ + uifi - uqUq , (3.17a) 
(p^ - a2P + ai)uo + ip- a2)uo,i + no,2 

= (p^ - a2P + ai)uo - (p- 02)^1,0 + U2,o + (p - a2)uoUo - iio?ti,o + uq^iUq , (3.17b) 
(g^ - a2q + ai)uo + (q - a2)iio,i + 1*0,2 

= (q"^ - OL2q + ai)uQ - (q- 02)^1,0 + ^2,0 + (q - a2)uQUQ - UQUifi + uq^iUq . (3.17c) 

Let us focus on ()3.17bp and (j3.17cp . By substraction one can delete uo,2) and to delete M2,o from 
the remains one can make use of (j3.7ap where we take ^ = 1, j = 0. Then, after some algebra 
we can reach to 

— G(—p,—q) G(—P,—q) . , :S > ^ 

— = \-{p + q- a2){uo - no) - uquq + uifi + no,i. (3.18) 

p-q + UQ-uo q-p 

This equation together with (I3.17ap composes of a closed-form system which is viewed as the 
extended three-component lattice BSQ equation. (We make a precise connection with the cor- 
responding system of [TT] in the next section). 
Next, by the transformation (if a ^ b) 

Sa,b = Sa,b - 7-— (3-19) 
— a 

we rewrite the relation (|3.8p and its hat-*; counterpart in the form 

(p - a)Sa,b - {P- b)Sa,b = VaWb , {q - a)Sa,b - {q- b)Sa,b = VaWb , (3.20a) 

PbSa,b - PaSa,b = SaWb - vjb + {p - a2)VaWb , (3.20b) 

qbSa,b - qaSa,b = SaWb " Vatb + {q - 02)^0 "Sfe . (3.20c) 

Then, from (]3.20bp ^ deleting % by using p.l4ap and by using the hat-g version of (|3.9ap and 
also making use of (|3.14bp we get 

[q - a)VaWb = 7 rr^ -, rr^VaWb + PaSa,b - PbSa,b ■ (3.21) 

[p - b)wb - (g - b)wb 
Using ()3.20ap once again we can rewrite (I3.2ip to the following desired form 

^ :^WbVa-^WbVa G(-a,-b) ~ 

VaWb = Wb-r -rr^ 7 TT^ + 7 TT ^^a,b ■ (3.22) 

[p - b)Wb - [q - b)Wb [p-a){q-a) 

The system composed of p.22p and p.20ap can be viewed as the extended three-component 
lattice MBSQ/SBSQ system, in the sense that the extended versions of both the 9-point lattice 
MBSQ of [23] as well as of the 9-point Schwarzian BSQ lattice can be obtained from this system 
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by elimination of two of the three variables. An (equivalent) alternative form of the extended 
three-component lattice MBSQ/SBSQ equation consists of (j3.20ap together with the following 
equation 

f^,WbVa-^,W,Va G{-a, -b) ^ 
VaWb = Wb-, TT^ -, + 7 Tw 7T'5a,6 • (3.23) 

[p - b)Wb - [q - b)Wb [p-o)[q-b) 

(|3.23p can be derived similar to (j3.22p . We note that (j3.22p and (|3.23p can be transformed to 
each other in the view of (j3.20ap . In fact, noting that the relation G{—p, —b) — G{—p, —a) = 
G{—a, —b), a substraction of (|3.23p from (j3.22p gives 

^ ^ WbVgWb WiVaWb 

^a,b ^a,b {p—a){p—b) {q-a){q-b) 



{p-a){q-a) {p - b){q - b) {p - b)wb - (q - b)wb 
which, by replacing VaWb and VaWb using (I3.20ap . further reduces to 

{p - a)wbSa,b -{q- a)wbSa,b 



{p - b)wb - (q- b)wb 



and this is an identity as is evident from (|3.20ap . 

Finally, to derive the coupled system (|1.4p we rewrite (j3.2ip in the form 

/ s <^ PbVaWb - qbVaWb - , ? a 

[q - a)VaWb = ^ VaWb + PaSa,b " Pb^a,b 

{p - b)VaWb - (g - b)VaWb 

_ Pb[{p - a)Sa^b - {p- b)Sa^b] - qb[{q - a)Sa,b - (g - b)Sa^b] 
{P -b){q- a)Sa,b - {q- b){p - a)Sa,b 

X [{q - a)Sa,b -{q- b)Sa,b] + PaSa,b - PbSa,b , (3.24) 

and introduce the variable u by setting 

S ^ _( {p-a)Pb y( iq-a)Qb Y^^ .325. 



(3.26) 



Furthermore, (jl.Sp follows from (|1.4p by making use of the equality 

jVaWb) _ VgWb VgWb 

(7;„Sb)~" W I I ■ 



We have now all the ingredients in place to make the identifications with the systems exhibited 
in[n]. 



4. Deformation of the extended lattice BSQ-type equations 

In this section, we identify the extended closed-form lattice ESQ equations obtained in the 
previous section, with the four extended lattice BSQ-type equations (A-2), (B-2), (C-3) and 
(C-4) of ref. [11] by means of simple point transformations. 
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Case A-2 : 

Starting from (|3.15p and introducing point transformation 

X 1 

Va = , UQ = Z - ZQ, Sa = {v - VaVa) , (4.1) 

where 

Xa = {p-ar^iq-ar'^ci , (4.2a) 

zq = (c3 - p)n + (c3 - g)m + C2 , (4.2b) 

Va = Xa{zo - Cs) , (4.2c) 

and ci, C2, C3 are constants, (j3.15p yields 

y = zx — X, y = zx — X , (4.3a) 

~ , , -G{-p,-a)x + G{-q,-a)x 
y = xz — oqx H — , (4.3b) 

z — z 

where bo = —02 + 3c3. This is the (A-2) equation given in ref. [TT] and the parameter 60 can be 
removed by a transformation 

Alternatively, employing the point transformation 

X 1 

Wb = —, Uq = z- Zq, tb = —{y - Wbyb) , (4.4) 

Xb Xb 

where 

Xb = i-p + bn-q + brci , (4.5a) 

zo = -(c3 + p)n - (c3 + q)m + C2 , (4.5b) 

Vb = Xb{zo - C3) , (4.5c) 

and ci, C2, C3 are constants, we find 

y = zx — X, y = 'zx — X , (4.6a) 
V = ^z- 6'o§ + -Gi-P,-b)x + G{-q,-b)x 

where 6q = a2 + 3c3 . This equation is related to (|4.3p by the transformation 

p — )■ —p, q —q, n — > — n, m —m, 02 — >■ —02, b —a, 
which corresponds to the reversal symmetry of (A-2), (see |llj). 
Case B-2 : 

Now we look at the extended three-component lattice BSQ equation (|3.17ap and (|3.18p . Making 
transformation 

X = uq — xq , (4.7a) 
y = ui,o - xqUq + yo , (4.7b) 

Z = Uo,l - XqUq + Zq , (4.7c) 
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with 

xq = np + mq + ci , (4.8a) 
yo = ^{np + mq + cif + ^{np^ + mq'^ + C2) + C3 , (4.8b) 

zq = ^{np + mq + ci)'^ - ^{np'^ + mq'^ + C2) - C3 , (4.8c) 

and constants Cj,{j = 1, 2, 3), we have 

z = XX — y, z = XX — y , (4.9a) 

~ ~ f~ \ G{-p,-q) 
z = XX — y — a2\x — x) — ai — — — — . (4.9b) 

X — X 

This is the (B-2) in Ref.|llj and the parameter ai can be removed by transformation 
Case C-3 : 

Next, we come to the extended three-component lattice MBSQ/SBSQ equation (j3.20ap and 
(|3.22p . By the transformation 

p — a\ I q — a 



Va = (p-ariq-ary, (4.10b) 

Wb = ip-b)-''{q-b)-"'z , (4.10c) 

then we have 

X — X = yz , X — X = yz , (4.11a) 

~ -G{-p,-b)zy + G{-q,-b)zy - u 

yz = z — hG[—a,—b)x. (4.11b) 

z — z 

Similarly, by the same transformation from (j3.20ap and ()3.23p we have 

X — X = yz , X — X = yz , (4.12a) 

~ -G{-p,-a)zy + G{-q,-a)zy 10K^ 
yz = z — h Ct(— a, — ojx . (4.12b) 

z — z 

This gives (C-3) equation in [11]. ()4.12p and (j4.1ip are related by 

n —n, m — )■ —m, y ^ z, z ^ —y, a ^ b 

which gives the reversal symmetry of (C-3) [11]. Besides, as we mentioned before, (|4.11bp and 
(|4.12bp are also connected through 

zy-zy - 

Z^z — = X — X, 

z — z 

which holds in light of (j4.11ap . 
Case C-4 : 

The identification of Case (C-4) of [Tl] is somewhat indirect and relies on the important obser- 
vation that ()4.1ip and (j4.12p share the same solution (j4.10p . Thus, from these realtions we have 
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where 



— X = 


yz , X — X 


= yz, 




yz = 


Pa,b zy - 
z — — 

z — 


Qa,bzy 

Z 


,b{x + x) 


Pa,b = 


-1{G{-P, 


-b) + G{-p, 


-a)), 


Qa,b = 


-l{G{-q, 


-b) + Gi-q, 


-a)), 


Ga,b = 


G{-a,-b) 
2 







Consider the following transformation 



xi - Ga,b yi 



^y= ,2= tl . (4.15) 

'^Ga,b{xi+Ga,b)' Xi+Ga^ Xi + Ga,b 

Imposing this transformation on ()4.13p . ()4.13ap yields 

xi-xi = yizi , xi-xi = yizi, (4.16) 

which further gives 

- XlZl - XiZl 

xi = — . (4.17) 

zi - Zi 

Meanwhile, ()4.13bp multiplied by {Ga,b + xi){Ga,b + xi) yields 

yiti = A + xJi - 6, (4.18) 

where ^ 

^ _ {GgM + Xi)zi{Pa^bZiyi - Qa^bZiyi) 
{xiZl - XlZl) + Ga,b{zi -Zl)' 

Substituting (j4.17p in to A we then have 

~ Pa,bZiyi-Qa,bZiyi - 2 on^ 
yizi = Zl ~ _^ h xixi - (4-20) 

Zl Zl 

and ()i:20D constitute the (C-4) equation of [TT] . 

In conclusion, we have in this section identified all deformed equations of extended lattice 
BSQ type of |llj with equations in section 3 arising there from the direct linearization scheme 
of section 2. We note that, whilst the results in [11] were obtained from a systematic search 
programme, the equations that are found from the classification conditions come somewhat as 
separate cases through the elimination of choices and the implementation of certain symmetries. 
However, search does not really reveal the interconnections between the variables of the different 
cases, and, in contrast, that is what the DL approach provides. Thus, deriving all equations from 
the single framework all the interconnections between the variables of the various systems are 
revealed. In the next section we show how explicit soliton solutions are obtained systematically 
from the DL scheme outlined in section 2 and applied to the choices of variables associated with 
the extended BSQ systems. 



(4.13a) 
(4.13b) 

(4.14a) 
(4.14b) 
(4.14c) 
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5. SOLITON SOLUTIONS 

In the case of soliton solutions, we start from the fohowing exphcit expression for the infinite 
matrix C 

3 

j=l j'=l 

which can be obtained by choosing in ()2.ip the integration contours Tj in the complex plane 
of the variable k to be Jordan curves surrounding simple poles at locations kjj/ in the complex 
plane, with residues Ajj/ of the measures d\j{k), assuming the latter to be meromorphic in k. 
In this case we can derive the following expression for the infinite matrix U 

3 Nj 

i=i j'=i 

in which the infinite-component vector obeys the linear equation: 

3 ^^7 

+ L L k-u,{k,,) = p'""' ■ (^-^^ 

Eqs. ()7.7p and ()5.3p are derived from the basic definitions in the Appendix. 

By setting in eq. (I5.3P k = kj/^ii, where j' = 1, 2, 3, i' = 1, . . . , Nj', the equation becomes a 
linear system for the quantities Uk. ., , with a finite block-Cauchy matrix of the form: 

which is a 3x3 block matrix with rectangular blocks of size NiXNj. We use the capital compound 
indices I = and J = (j, j') to simplify the notation. 

We can make the soliton solutions now explicit by assuming that the coefficients j/ are 
chosen such that the matrix 1 + M is invertible, in which case we can make the following 
identifications for the quantities introduced in (j2.9p and (j3.6p respectively: 



+ M)-^ K' r , Sa,b = {-bl + K')-^ {1 + M)-^ {al + K)-^ r (5.5) 
in which K and K' are the block diagonal matrices: 

K = diag{ki^i,...,ki^Ni;k2,i,...,k2,N2;k3^i,...,k3^N3) , (5.6a) 
K' = diag{-uji{ki^i),. . . ,-uJi{ki^Ni);-^^2ik2,i),- ■ ■ ,-i^2{k2,N2);-^3ik3,i),- ■ ■ ,-^^3{k3,N3)) , 

(5.6b) 

and where r and *s are the vectors 

r = diag(^pki^^,...,pkij^^;pk2,i,---,Pk2,N2'^Pk3,i,---,Pks,N3) > (5-7a) 

*s = diag j), . . . jyj;a_^2(fc2,i), • • • ,f^-a;2(fc2,Ar2);'^-'^3(fc3,i) • 

(5.7b) 
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The soliton solutions of the other relevant variables, defined in (I3.5p are obtained in a similar 
way, and lead to the expressions 

Va := 1 - 's {1 + M)~^ (al + K)-^ r , Wb := 1 + 's {-bl + K'yH^ + M)"^ r , 

Sa := a - *s K' {1 + M)~^ {al + K)-^ r , := -6 + *s (-61 + K')"^ (1 + M)"^ K r , 

ra := - K'^ {1 + M)-^ {al + r , Zb := + {-bl + K')~^ {1 + M)"^ K'^ r 

Note that all these expressions (15.5P and (jS.Sp can be made quite explicit after choosing the 
various parameters of the solution. Finally the relevant r-function associated with these soliton 
solutions is given by 

r := det (1 + M) (5.9) 
and it can be shown that it obeys the following trilinear equation: 

(3p^ — 2a2P + ai)TTT + (3q^ — 2a2q + ai)TTT — {p — q)'^TTT 

= [{p^ + pq + q^) - a2{p + q) + ai](^TT+jjT) . (5.10) 

This trilinear equation is derived from a combination of (j3.13p . (j3.15p and (|3.12p together with 
the identifications 

T T T T T 

Vp=- , Wp=-, Vg = -, Wg = -, {q-p)Sp^g=-. (S.ll) 

T T T T T 

The explicit expressions (|5.8p can be subsequently inserted in the formulae of the previous 
section identifying the various dependent variables of the systems of [H], and thus we get the 
explicit multi-soliton solutions for the latter. Furthermore, the identifications (|5.1ip can be used 
to express these quantities in terms of the r-function. It is to be expected that the trilinear 
equation could be cast in bilinear form, but at the expense of having to introduce additional 
r-functions (e.g. by introducing functions obtained by applying a shift on (|5.9p in a lattice 
direction other than the ones appearing in the equation itself), but we will not pursue that line 
of investigation here. 



6. Lax representations 

The relations (13. 4p are the basis for the derivation of Lax pairs for the lattice equations list 
in section 3. In fact, by selecting specific components of these vectors, or combinations thereof, 
we can constitute the basic vector functions in terms of which we obtain the relevant linear 
problems. This the DL scheme will provide some of the Lax pairs directly. In some cases we 
find it more convenient to rely on the 3D-consistency of the systems to produce the Lax pairs, 
along the lines of the papers [Mill], cf. also (30l[28] for the BSQ case. 

Case A-2 : 

For the (A-2) equation (I3.15p . we define an eigenvector of the form(cf. |24j ) 



0=(((a + A) iufc)o, (ufc)o, (Aufc)o) 
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Here for a infinite-component column vector v, (v)o stands for the 0-th component of v. In that 
case one can derive the fohowing Lax relation cf) = L^'^cf) with: 

Va 

L^^ = \ p-uo 1 1 (6.1) 




* p-a2 + ^ 

in which * = {p — uo){p — Q2 + Sa/va) — PaVa/va, and M-^"^ is obtained from (j6.ip by replacing 
p hy q and ~ by ^. The compatibility condition of the Lax pair, i.e. the matrix equation 
j^A2j^A2 ^ j^A2j^A2 y^^^^ ^j^^ equations (|3.13ap and (|3.13bp . together with the relation 

p- = {p- a)^ - {q- a)^ , (6.2) 

Va Va Va Va 

The latter relation follow readily from (j3.9ap . 
Case B-2 : 

For the extended lattice ESQ equation (j3.17ap and ()3.18p . then one can define vector 

ip = {{uk)o, (Aiifc)o, (wi(-A)w2(-A)'Ufc)o)'^, 
then from ()3.4p . we obtain 

^ = L^SQ(f) , (6.3a) 
^ = M^^^ct) , (6.3b) 

where 

/ p-uo 1 \ 

lBSQ ^ _~ _ p^^^ 1 I (g_4) 

\ * -1*0,1 + 202-^0 - p - 2a2 + uq ) 

in which * = G(/c, -p) - {p - uo)[{p - a2){p + uo) + no,i] - {p + uq- 2a2)(ni_o + ai), and where 

jyjBSQ jg obtained from (|6.4p by replacing phj q and ~ by ^. The compatibility L^^Q j^BSQ _ 

jyjBSQj^BSQ Y^g^^s to the equations 

Si,o - Si,o = (p-q + UQ- uo)uo - puo + quo (6.5a) 

and 

uo,i - So,i = {p - q + uo - uq)uo - puo + quo, (6.5b) 

together with eq. (j3.18p . These three equations constitute a system which is equivalent to the 
extended lattice ESQ equation (I3.17ap and p.lSp . 

In the following cases, we give Lax pairs which are derived on the basis of the 3D consistency 
property for the systems of three coupled equation which we derived in section 4. In the generic 
case for a 3-component system containing variables (x, y, z), we would need to set {x = g/ f, y = 
h/f, 1 = k/ f) and constitute the vector is 4> = {g, h, k, f)'^ leading to a Lax pair in terms of 
4x4 matrices (see, e.g., (j6.23p and (j6.3ip below). However in some cases the 4x4 situation 
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simplifies to a Lax pair in terms of 3 x 3 matrices for the 3-component system (e.g., in the cases 
of the A-2 and B-2 systems). 

Case A-2 revisited: 

For the extended three-component lattice equation (j4.3p . we extend the lattice equation into a 
third dimension by introducing a new variable / associated with a new shift ~ and a new complex 
parameter r 



y = zx — X , y = zx — X , y = zx — x , 

-G{—p, —a) X + G{—q, —a) x 



y = xz — box + 



y = xz — bQX + 



y = xz — Bqx + 
From ()6.6ap . it is easy to derive 



z — z 

-G{—p, —a) X + G{—r, —a) x 
'z — z 

-G{—r, —a) X + G{—q, —a) x 



X 



x — X ~ X — X 
, X 



X 



z — z 
xz — xz 
1, — z 



z — z 
xyz — "xz 



X — X 
z — z 

^ xz — xz 



y 



z — z 



(6.6a) 
(6.6b) 

(6.6c) 

(6.6d) 

(6.7a) 
(6.7b) 



By the 3D-consistency, setting x = z = y = then the Lax pair is 



011 



/ 



= -Lii^ii , = Mu4>u > 
I z 



'11 



z 

X 



in which 

-1 0\ 

-1 

(|+5„)j_G(zl^ -(| + 6o) 

V -F X / 

where the matrix iWii is the hat-g version of Li\\. The compatibility gives the relations: 

^ X — X xy — xy 



(6.9) 



(6.10) 



z — Z X — X 

together with eq. (I4.3b|) . The two equations (I6.10p are consequences of the 3-component system 
which constitutes the case (A-2). We note that there is a zero column appearing in Ln as well 
as in iWii, which indicates y is not necessary to be used to generate a Lax pair. Let us remove 
the last columns and rows in Ln and iVfn, and examine the remains: 



'12 



'12 



in which 





( ^0 ^ 




012 = 








{^2 J 





-/W"l20l2i 





z 

X -1 

I br,)z G{-p~a)x G{~r~a) 



-1 





(6.11) 



(6.12) 
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where the matrix -M12 is the hat-g version of L12. The compatibihty of ()6.1ip gives ()6.10p and 
(f43bl) as well. 

Case B-2 revisited: 

Similarly, for the extended three-component lattice BSQ equation ()4.9p . by considering the third 
dimension of this lattice equation, then we have 



Z = XX 


— y, = XX — y 


z = xx - y , 


(D.loaj 


Z = XX 


-y - OL2{x - x) 


Oil _ _ , 
X — X 


(6.13b) 


Z = XX 


- y - a2{x - x) 


G{-p,-r) 
ai _ ~ , 

X — X 


(6.13c) 


Z = XX 


-y - a2{x - x) 


G{-r,-q) 
— cxi — . 

X — X 


(6.13d) 


From (|6.13ap. it is easv to derive 








^ y- 

X = — — 
X — 


y ^ y-y 

^ 1 X ^ — , 
X X — X 


- y-y 

x = _ ^ , 
z — X 


(6.14a) 


- xy- 
z = 

X - 


-xy ^ xy — 

^ , z= _ 
- X x — 


xy ~ xy — xy 
X X — X 


(6.14b) 


By the SD-consistency, setting x = 


00 ' ^ ~ 00 ' ^ ~ 


then the Lax pair is 




4>2l 


= -^^21021 ' $21 


= -^"21021 , 


(6.15) 



0: 



21 



21 



(6.16) 



in which 

fx -1 \ 

y -10 

* z — a2X + ai a2 — X 

\ y -X J 

with * = —{z — a2X + ai)x + y{x — Q2) + G{—p, — r), and where the matrix -M21 is the hat-g 
version of L21. The compatibility yields the relations 

- y-y z-z 



f <Po\ 

02 



X = - — - 6.17 

X — X X — X 

together with the relation ()4.9bp . Similar to the case (A-2) the Lax matrices L21 and M21 can 
degenerate to 3 x 3 case by removing their last rows and columns. The results are 



022 = -^22022 1 022 = -^22022 ; 

-1 









22 = 




, L22={ 




{ <P2 ) 





(6.18) 
(6.19) 



in which 

X 

y -1 

* Z — a2X + Qfl 02 — X 

with * is defined as in (j6.16p and the matrix iW22 is the hat-f; version of ^22- The compatibility 
yields the relations (|6.17p and ()4.9bp 
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Case C-3 : 

For this case, which includes the extended three-component lattice MBSQ/SBSQ equation 
(|4.1ip . we are led by considering a third dimension of this lattice equation to the following 
relations: 



X — X = yz , X — X = yz , x — x = yz 
-G{-p,-b)zy + G{-q,-b)zy 



yz = z- 
yz = z- 
yz = z- 



z — z 

-G{-p, -b)zy + G(-r, -b)zy 
z — z 

-G{-r,-b)zy + G{-q,-b)zy 



z — z 



+ G{—a, —b)x 
+ G{-a, -b)x 



+ G{-a, -b)x. 



From ()6.20ap . it is easy to derive 



^ xz — xz — 
X = — —— , X 



z — z 

X — X 



XZ — xz ^ 
— — — , X 

z — z 



xz — xz 



z — z 



y 



y — y — X — X 
-z- — - , y = - — - 



z — z z — z z — z z — z 



z — z 

_ ^2 



y — y ^ X — X 
-z- — = , y = - — - 



z — z 



y-y 



z — z 



By the SD-consistency, setting x = ^, y = ^, z - ^ 
in which -0]^ = (-00, i>i, i>2, i's)^ and 






— then we have the Lax pair 



1 

X 





where the matrix is the hat-g version of -L3. The compatibility yields now: 





( 


—z 





1 







—z 


z 




zy 












-G{-a, 



(6.20a) 
(6.20b) 

(6.20c) 

(6.20d) 

(6.21a) 
(6.21b) 

(6.22) 
(6.23) 



xz — xz 



X 



y 



-Z-, 



y-y 



(6.24) 

z — z " z — z 

together with the relation ()4.11bp . They constitute the case (C-3) of the extended BSQ family. 

Alternate Case C-3 : 

In a similar way as in the previous case, for the extended three-component lattice MBSQ/SBSQ 
equation ()4.12p . by considering a third dimension where we impose this lattice system, we have 

X — X = yz , X — X = yz , x — x = yz , (6.25a) 
-G{-p, -a)zy + G{-q, -a)zy 



yz = z- 
yz = z- 
yz = z- 



z — z 

-G{-p, -a)zy + G{-r, -a)zy 
z — z 

-G(-r, -a)zy + G{-q, -a)zy 



+ G{—a, —h)x , 
+ G{—a, —b)x , 
+ G(-a, -b)x . 



(6.25b) 
(6.25c) 
(6.25d) 
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Obviously, the identities (I6.2ip are still hold. By using SD-consistency, and setting x = y 

Mo ' 



l^, then from ([OT|) . (|6.29c|) and ()6.29dp we have the Lax pair 



in which i/'2 = (/^o, A^i, Ai2, /"s)^ and 





/ 


—z 








1 


\ 


1 







—z 









^; 




zy 





—z 









v 


-Gi-a,-b)f 





-G(-r,-a)f 


G(-p,-a)f + G(-a, 





(6.26) 



(6.27) 



where the matrix M4 is the hat-q version of L4. The compatibility now produces the relations 



xz — xz 



xy — xy 



,'y = -zy^ (6.28) 

z — z y — y z — z 

together with ()4.12bp . Once again this system is a consequence of the original lattice system. 
Case C-4 : 

For this case, which includes the extended three-component lattice MBSQ/SBSQ equation (14.16P 
and (|4.20p . we have the following SD-consistent relations: 



Xl - 


Xl = 


yizi , Xl - 


-xi= yizi , Xl - 


- Xl = yizi , 


(6.29a) 


yizi 


= Zl 


Pa,b Ziyi - 


Qa,bZiyi , - 
h XlXi 


~ '^a,b ; 


(6.29b) 






Zl - 


Zl 


yizi 


= Zl 


Pa,bZiyi - 


Ra,bZiyi , 
h XlXl 


^2 


(6.29c) 






Zl - 


Zl 




yizi 


= Zl 


Ra,b ZlVi - 


QafiZiyi , - 
— 1- XiXl 


^2 


(6.29d) 



where Ra,b = -i(G(-r, -5) + G(-r, -a)). From (iOgil) we have same consistent forms as (|6.2ip 



with {x,y,z) replaced by {xi,yi, zi). Setting xi = 
in which rp^ — i^Oi ^i^ ^2, 1^3)^ and 



1 

Zl 



I/O' ^1 



M5-0 



3' 



we have the Lax pair 
(6.30) 





-Zl 








1 







-Zl 





Xl 




ziyi 





-Zl 





v 


'^a,byi 


XlZl 

yi 


rLa,b 


r-a,b y^ 



(6.31) 



where the matrix is the hat-g version of -L5. Now the compatibility yields 



Xl 



xiyi - xiyi 



Xl 



XiZi — XiZi 



yi 



-Zi- 



yi - yi 



(6.32) 

2/1-2/1 Zl- Zl ' - Zl- Zl 

together with the relation (|4.20p . They constitute the case (C-4) of the extended BSQ family. 
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7. Conclusions 



Since the realization that the notion of multidimensional consistency is a key criterion for 
integrability of partial difference equations, the search for new integrable equations on the 2D 
lattice using this notion as the defining property has led to many novel results. The ABS 
classification, (2], of single- variable quadrilateral lattice equations has demonstrated that such 
systems are highly restricted (up to equivalence), but the non-scalar case remains still to be 
largely explored. The search in of 2-and 3-component systems of a form inspired by earlier 
results from the BSQ family of lattice equations, [24^ [T9l |28] seems to confirm the assertion 
that the realm of such discrete equations has only few members, but the examples are rich and 
important. The parameter-extensions of the BSQ systems established in are noteworthy 
because the extra parameter freedom allows their soliton solutions to exhibit some new features, 
cf. |12| I13j. On the other hand, the "search methodology" exploited in [Tl] does not reveal all 
there is to know about these equations, and this is supplemented by the structural approach 
wich we took in the present paper. Thus, by extending the DL approach of [23] to include more 
general types of dispersion relations, we have been able in the present work to identify all the 
variables and their interconnections of the extended BSQ systems of [11]. In fact, this amounts to 
establishing Miura type (i.e., non-auto Backlund) transformations between these various systems 
(A-2), (B-3), (C-3) and (C-4), which involve also point transformations. Furthermore, the DL 
structure can be used effectively to yield explicit multi-soliton type solutions for all these systems 
by specifying in an appropriate way the contours and integration measures which define some 
of the main objects (e.g., in (j2.ip ). For more general measures and contours the DL also can 
yield wider classes of solutions, such as inverse scattering type solutions or scaling reductions. 
To explore further, and make explicit, the latter line of work on the basis of the present results, 
will be the subject of future investigations. 
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In this Appendix we derive some integral formulae from the general abstract structure dis- 
played in section 2. Starting from (j2.ip . i.e.. 
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is equivalent to the defining relation 

U = C-UftC. 

In fact, 



N 

7=1 •''^1 



Using (|L2|) it is then 



N 



= -unc. 

On the other hand, starting from (j7.ip and using (|7.2p . one has 

C = / ^^i(^) c*_^^,(;^)CJ_^^(fc) +Un^ dXj{k) pk ck c*_^^.(fc)CJ- 

N 

= 2^ / ^^i(^) '"fc c*„^^.(;.)(j„^^(fe) + ?7 n C, 

which yields ([73]) if (IZ3I) holds. 

Then, quite similar to the derivation in (j7.5p . if starting from (jS.ip . i.e., 

3 

and using ()7.2p and ()7.4p . we immediately have 

3 

j=i j'=i 

which is (j7.7p . Finally, for ti^ we have 

\ ^ \ ^ , + K — tOj[Kj ji) 

3 -^^j 



which gives ()5.3p . where we have successfully made use of 

A Cfc = A; Cfc, 4, *A = /c' 4,, *An + OA = O, O = e *e 

and 4' O Cfc = 1. 
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